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Cooper Pair Transport in a Resistor-Biased 
Josephson Junction Array 

Sergey V. Lotkhov, Vladimir A. Krupenin, and Alexander B. Zorin 

Abstract 

The dc transport properties of long arrays of small Al Josephson junctions, biased through on-chip 
Cr resistors, are studied. The /^-characteristics show a large Coulomb threshold for current as well 
as negative-slope regions indicating the regime of autonomous Bloch oscillations up to rather high 
frequencies of / = J/2e ~ 1 GHz, comparable to those reported by other groups for single junctions. 
On the other hand, a small depth of the back-bending implies a low duty cycle and a broad spectrum 
of the oscillations, which we attribute to the insufficiently high impedance of the bias resistors. A 
self-sustained switching process at a small bias current is used to study the statistics of the switching 
voltages and to determine the effective Bloch capacitance which was found to considerably exceed the 
geometric junction capacitance. 

Index Terms 

Charge transfer, current, Josephson arrays, SQUIDs, superconductor-insulator-superconductor de- 
vices, thin-film devices, stripline components. 

I. Introduction 

The application of single-charge tunneling effects in electrical metrology has been considered 
in a view of capacitance and current standards based on charge quantization phenomena (see, 
e.g., review [1]). The proper charge quantization in circuits with small tunnel junctions has been 
achieved for low junction transparencies: R T 3> Rq = h/e 2 « 25.8 kf2 for the single-electron 
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tunnelling resistance, and Ej <C E c for Cooper pair devices, with Ej oc R^ 1 and E c being the 
Josephson coupling and charging energies, respectively. For example, remarkable advances have 
been made in development towards the single-electron standard of capacitance [2]. The accurate 
operation of these circuits, based on pumping of single electrons, has been, however, limited 
to small currents, I ~ 1 to 10 pA, because of the stochastic nature of single electron tunneling 
with a relatively long time constant. 

Here we address the Cooper pair transport in the case of substantial Josephson coupling Ej > 
E c , i.e., when both the Josephson phase and the charge exhibit significant quantum uncertainties. 
In such rather transparent, and therefore large-current-capable junctions, high accuracy could still 
be achieved due to the exact 2e-periodicity of their energy bands as a function of the quasicharge 
[3]. 

II. Basic idea and sample 

The related transport phenomenon is known as Bloch oscillations in a current-biased junction 
with a fundamental current-frequency relation I = 2ef [3]. Successful attempts to phase-lock 
the Bloch oscillation by an external signal were made in the early 1990s (see, e.g., [4], [5], 
[6]) in the single junctions biased through high-ohmic on-chip microresistors for frequencies 
up to / ~ 10 GHz. Unfortunately, the linewidth of oscillations was found to be rather large, 
Sf ~ 1 GHz, which was supposedly due to significant thermal fluctuations which smeared the 
current plateaus. Moreover, a possible size of these plateaus was expected to be rather small, 
scaling as a small threshold voltage in a single junction. 

Our approach is based on ID arrays of small Josephson junctions. As was shown theoretically 
for current-biased arrays [7] and proved experimentally for very long voltage-biased Al arrays 
(N ~ 200) [8], these arrays are analogous to long Josephson junctions: Similar to the fluxons 
formed in the long junctions, the charge profile in the array is governed by the sine-Gordon type 
of equation, describing the Cooper pair solitons of size (in number of junctions) M ~ \Jc/C < 
N [8]. Here we denote as C the effective junction capacitance, related to the curvature of the 
ground state energy E(q) and, hence, the shape of this 2e-soliton; C <C C is the self-capacitance 
of the intermediate islands. The coherent motion of a train of 2e-solitons along an array biased 
via high-ohmic resistors can be described by a set of non- stationary equations: 
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L B (q i )^ + ^^-=^-(q i ^-2q i + q i+1 ), i = l,2,... (1) 

which includes the quasicharge-dependent Bloch inductances of the junctions Lb(<Ji) which were 
recently introduced in [9], with corresponding initial and boundary conditions. Here we denote 
the set of junction quasicharges as {q^t)}. In the limit of small Josephson coupling, A = 
Ej/Eq <C 1, the periodic term dE(qi)/dqi takes the form of the sawtooth function with the 
amplitude e/C, where C is the capacitance of the individual junctions in the array. In the 
opposite case, A ^> 1, there is an analytic expression [3]: 

dE{qj) e . 

dq = a(A)— sm(7rgi/e), (2) 

where a(A) = 2 11//4 7r 1 / 2 A 3 / 4 exp — (8A) 1 ' 2 is a numerical factor describing exponential sup- 
pression of the Bloch band width for A — > oo. 

In this paper, we report the dc properties of autonomous arrays of iV = 60 SQUIDs (see 
the SEM photo in Fig. [l]), fabricated using the Al/oxide/Al-Cr shadow evaporation process (see, 
e.g., [10]) in a 4-point layout with Cr microstrips of 50 /im-length with R Cv w 550 kf2 each. The 
tunnel resistance was ~ 8kf2 and an estimated capacitance was Ct ~ 0.36 fF per junction, 
which corresponds to the maximum Josephson coupling energy of each SQUID, Ej rs 160 /zeV 
and the charging energy E c = e 2 /4C T ~ 110/ieV per link, yielding their ratio A w 1.5. We 
roughly estimate C ~ 50 aF. 

The samples were measured in the current-bias mode with an external high-ohmic bias resistor, 
i?B = 100 MQ. The capacitance of each connection line to ground was Cl ~ 1 nF, resulting in 
an i?C-constant of the bias circuitry of td ~ RbCl = 0.1 s. Our microwave-tight holder was 
supplied with Thermocoax® filters (jacket 0.35 mm, 110 cm long), anchored to the mixing 
chamber. 

III. Results 

Several important peculiarities can be observed in the /^-characteristics of the arrays shown 
in Fig. |2| In our opinion, the most remarkable transport property of the system relates to the 
bias current range Jbi as < 300 to 350 pA, with zero to negative slope of the IV-curve (see 
the blow-up of the corresponding region in the right inset in Fig. 0), which we assume to 
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be an indication of the Bloch oscillation regime [3]. Due to the moderate current values, the 
ramping of quasicharge in this transport regime occurs sufficiently slowly that the system always 
stays in the zero Bloch band and experiences 2e-periodic oscillations of the voltage with the 
frequencies up to / = J bias /2e ~ 1 GHz. At larger currents, the average voltage across the array 
gradually increases, presumably due to excitation of the upper energy bands [3], thus setting a 
high-frequency limit for possible application of these particular Josephson arrays as sources of 
a quantized current. Note that, despite a large number of junctions in our arrays, this limitation, 
being at the level of GHz, is virtually of the same scale with the frequencies of phase-locking 
experiments on single junctions [4], whereas the typical voltage scale is much larger, which is 
promising for a better observability of the Bloch oscillations. 

On the other hand, small depth of the back-bending, typically, few percent, points out the low 
duty cycle of Bloch oscillations with the short pulse duration of r arr < Rc T C &rr ~ (5 GHz) -1 <C 
f~ x , where we use a rough estimate for the "input" capacitance of the array, C arr ~ {2CCq) 1 ^ 2 ~ 
0.4 fF (see below for our estimation of the effective value C). We attribute the present form 
of oscillations to insufficiently high impedance of the biasing resistors. Simple estimations 
show that the observed upper frequency limit for a single-band behaviour of about 1 GHz is 
achieved by the system biased only slightly above its Coulomb blockade threshold, therefore 
resulting in a typical decay time profile of oscillations. A broad spectrum of such oscillations, 
with considerable contribution of higher harmonics, makes it difficult to ensure their effective 
phase-locking to an external high-frequency signal. In particular, the preliminary experiments 
on microwave irradiation of our resistor-biased arrays up to frequencies of about 1 GHz did not 
result in observable features in the iV-curves. On the other hand, it was shown in [11] for 
the Josephson oscillations that even such oscillations may be effectively phase-locked using a 
driving signal of an appropriate waveform. 

As an indication of a hysteretical behaviour due to the inductance term in the equation 
of motion dU [8], [9], the biasing point at small currents J bias ~ 10 pA was found to be 
unstable, exhibiting relatively slow irregular voltage oscillations, see left inset in Fig. El As 
shown schematically in Fig. |3j random switching of the array from the blockade to the finite- 
current state was followed by a relatively slow recharging of the line capacitance Cl, which 
made possible the real-time observation of the switching cycles. In the blockade state, A^B, 
with a life-time At ~ r D , the current source charged the input terminal of the array until its 
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switching, B^C, to the finite-current state, /bias ~ 1 nA, followed by a rapid discharging, C^D, 
down to voltages V r ~ 500 //V and retrapping, D— >A, facilitated by the presence of fluctuations 
(cf. the thermally enhanced retrapping process in a shunted Josephson junction [12]). 

Whilst at T < 100 mK the average escape time, £ mcan = (At), was nearly temperature- 
independent (cf. [13]), as shown in Fig. |4j the switchings at T > 100 mK were due to thermal 
escape over the barrier whose height we have evaluated to be AU ~ 70 /JeV, as estimated from 
the slope of the Ahrrenius plot, see inset to Fig. HJ 

Using the approach [14] which is valid, strictly speaking, only in the weak-coupling case 
A <C 1, we obtained a rough estimate of the effective junction capacitance and the length of a 
soliton in the array. It is possible to express the barrier height AU(V) through the threshold 
voltage Vt of soliton motion: 



Taking advantage of the narrow range of switching voltages at T > 100 mK, V ~ V T , we 
estimated V T « 830 /iV, a soliton energy E = AU(0) = eV T /2 = e 2 /(2CC Q ) 1 ' 2 « 420 /xeV > 
Eq, the effective Bloch capacitance C ~ 1.6 fF ^> Ct, and the size of a soliton M ~ 8. One can 
see that the effective Bloch capacitance greatly exceeds the geometric junction capacitance which 
is obviously due to the strong suppression of the Bloch band width, AE = E max (q) — E mm (q), 
for the values of A > 1 [see, e.g., the limiting case, A>1, described by equation ©]. The shape 
of a static 2e-soliton in the infinitely long array is shown in Fig. 13 as calculated by numerically 
solving the equation ([Hi in the stationary case dqi/dt = 0, for several representative values of 
the ratio Ej/E c . 



Our dc measurements pointed out the presence of autonomous Bloch oscillations in a resis- 
tively biased array of small Josephson junctions. Compared to the experiments with single Bloch 
junctions reported by other groups, similar Bloch frequencies are expected, whereas the Coulomb 
voltage threshold and, as a consequence, the amplitudes of oscillations are considerably larger, 
making a potential advantage of using the arrays for phase-locking experiments. However, to 
enable an external synchronization of these oscillations one should realize at least several times 




(3) 



IV. Conclusions and outlook 
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higher bias impedance, possibly in form of very resistive microstrips. A certain increase in the 
Bloch frequency could be expected for the junctions with a larger Josephson energy. But since 
the voltage scale is also affected by the values of Ej and E c , a detailed optimization is still 
important. Using the statistics of the spontaneous state switchings, we evaluated the effective 
junction capacitance and the length of a Cooper pair soliton. These data can be helpful for future 
device development. For practical applications, the influence of the background charges in the 
islands of the array on the soliton dynamics should also be investigated. 
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FIGURES 

Fig.l. SEM picture of a fragment of the resistor-biased array of two-junction SQUIDs fabri- 
cated with the three-angle evaporation technique. 

Fig. 2. The iV-curves of two slightly different samples. Left inset: a typical time trace of the 
self-sustaining switching cycles observed at low bias currents. Right inset: a close-up look of 
the negative-slope segments in the iV-curves (shown by horizontal arrows) corresponding to 
single-band Bloch oscillations. 

Fig. 3. An equivalent electrical circuit (a) and a cycle diagram in the iV-plane (b), both 
explaining the switching/retrapping process. For the sake of simplicity, the on-chip Cr resistors 
and the four-point connection layout of the experiment are not shown. The Josephson array is 
schematically represented by the hatched box. In the presence of fluctuations (external noise, 
thermal fluctuations, etc.) there is no stable bias point along the load line, which results in the 
self-sustaining voltage relaxation oscillations observed in experiment. 

Fig. 4. Mean lifetimes of the blockade state. Inset: Ahrrenius plot at high temperatures. 

Fig. 5. Calculated shape of a 2e-soliton presented in a form of the island charges Qi = q i+ i — 
qi = C Vi, where is a quasicharge of an i-th junction and Vi is a voltage on the i-th island. 
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